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Abstract 

We derive the dimensional non-perturbative part of the QCD effective action for 
scalar and pseudoscalar meson fields by means of chiral and conformal bosonization. 
The related structural coupling constants L5 and Lg of the chiral lagrangian are esti- 
mated using general relations which are valid in a variety of chiral bosonization mod- 
els without explicit reference to model parameters. The asymptotics for large scalar 
fields in QCD is elaborated, and model-independent constraints on dimensional cou- 
pling constants of the effective meson lagrangian are evaluated. We determine also the 
interaction between scalar quarkonium and the gluon density and obtain the scalar 
glueball-quarkonium potential. 
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1. Introduction 



The low-energy effective action for light mesons |]l| describes their strong and electroweak 
interactions by means of several structural constants |]2[ which contain the information about 
the dynamics of quark and gluon interaction - Quantum Chromodynamics (QCD). These 
constants can be evaluated from the analysis of experimental data [|2|, Q , thereby giving a 
possibility to examine QCD at low energies. 

For the calculation of coupling constants of the chiral lagrangian direct bosonization 
methods of low-energy QCD @]-[|l3] and quark models |]ll|-[^ have been developed (see 
reviews ||T7| , 0). These methods give good and stable numerical estimates for certain di- 
mensionless chiral coefficients, namely for Li, L2, L^, L4, Lq, Lg, Liq 
case and for Liq, Ln, L19 



,[|T5[-|T9[inthe5f/F(3) 
in the Ua{1) extension ||2^ of the chiral lagrangian. Meantime 
the chiral coefficients L5, L7, Lg, L14, L15, Lig depending explicitly on QCD order parame- 
ters (such as the quark and gluon condensate p^) were not so far calculated in a model- 
independent way (see however large-iVc relations between them in p2| 



Thus, it is one of 



the goals of our paper, to find new constraints on those constants, in particular L5 and Lg 



from basic properties of the QCD vacuum [^ . We analyze the common structure of effec- 
tive lagrangians for pseudoscalar and scalar mesons as they are derived in general by the 
low-energy QCD bosonization procedure, and impose on such lagrangians the conditions to 
reproduce the behavior of the non-perturbative QCD vacuum energy in the large mass limit. 
These conditions together with a minimal number of physical inputs allow us to find the di- 
mensional chiral coefficients, in particular L5 and Lg. In Chiral Pertubation Theory (ChPT) 
the constant L5 enters the determination of the ratio of weak decay constants L^ : Fk '■ F^j 
of pseudoscalar mesons. The constant Lg characterizes the K^ — K^ meson mass difference 
as well as the current quark mass ratio {2ms ~ ""^d ~ ""^u) '■ ("^d """ ""^u) 0- 

In a general QCD bosonization approach the dynamical chiral symmetry breaking (DCSB) 
is implemented in order to derive the effective meson lagrangian. The DCSB is typically sim- 
ulated by means of a quark momentum cutoff A and a quark spectrum asymmetry M in 
QCD [^, 1^ or a dynamical quark mass Md in chiral quark models [|Tl|, |T^, |T8[ and quark 



models of Nambu-Jona-Lasinio type [15H|17|, [19|, |^|- The low-energy meson lagrangian 
is obtained from the bosonization of the quark determinant by application of the derivative 
expansion combined with the expansion in inverse powers of M or M^- It contains func- 
tionals of meson fields and external sources of different canonical dimension with coupling 
constants proportional to the DCSB order parameters. We find that the very structure of 
the above-mentioned functionals is universal in respect to different bosonization schemes 
and the model-dependent information is collected in structural constants. Just this observa- 
tion together with QCD vacuum energy constraints entails the definite predictions on chiral 



coefficients. As an interesting result, in addition we derive the scalar glueball-quarkonium 
potential directly from the QCD generating functional. 

The starting point of our analysis is the conventional QCD generating functional for 
colorless quark and gluon currents, 



Z{p,V,A,S,P)^ 
'DG'Dq'Dq exp 



/ 



d X 



y 



p] -AGl^f + m + ^S + l,P)q 



;i) 



Throughout, we will work in the Euclidean space. The covariant derivative is given by 
P = il'^id^, - iGlX^ + l^^ + 75^^) with tr 



A"A = 1/2 5 in color space. The external 
isoscalar and isovector colorless sources p(x), V^{x\ A^(a;), S{x\ P{x) serve for the building 
of Green functions of colorless currents as usual. The gauge fixing and ghost terms are not 
quoted explicitly but are understood to be taken into account in perturbative calculations. 
The derivation of the scalar-pseudoscalar meson lagrangian from QCD by direct bosoniza- 
tion 0-|T^ or by bosonization of a quark model [|n|]-[^ can be described schematically as 
follows, 

Z{p,V,A,S,P) = jVGVqVqVh5{h-{G%f)eM-S{q,q,G-p,V,A,S,P)) 

~ fvhVEVU exp{-Seff{h,J:,U; p,V,A,S,P) , (2) 

where the averaging over gluons G^ and quarks q, q is approximately replaced by the aver- 
aging over chiral fields U as well as over scalar glueball h and quarkonium E fields. Next 
the effective action is expanded in derivatives of meson fields and external sources. The 
model design consists in the choice of collective bosonic fields describing light mesons and of 
non-perturbative parameters specifying the DCSB (the type and magnitude of a momentum 
cutoff, the spectral asymmetry or the dynamical mass). 

The paper is organized as follows: In the next section we discuss the relations between 
the structure constants of the phenomenological low-energy meson lagrangian with those 
of model calculations. We will show that the values of dimensional structural constants, 
especially, of L^ and Lg are rather model-dependent. In section 3 we therefore extend the 
QCD bosonization approach to the scalar sector including all possible dimensional coupling 
constants. Comparing with the large scalar field asymptotics of the QCD vacuum energy, 
which is derived in section 4, we find additional constraints on the chiral structure constants 
which are discussed in section 5. In addition, this yields the scalar glueball-quarkonium 
potential. A summary of the results and a conclusion are given in section 6. 



2. QCD bosonization in the pseudoscalar sector 

In the chiral field sector {h, S ~ const) the effective lagrangian has the following structure, 

Seff = / d'^x (£2 + C4) + Swzw , (3) 

Here, the Weinberg lagrangian which is of chiral dimension 2 is given by 



C2 = -ft! [{D.UyO'^U - {x^U + U^x 



(4) 



where Fq ~ 90 MeV is the (bare) pion-decay constant, U is the usual S't/i? (3)-chiral field 
describing pseudoscalar mesons, UW = 1, det U = 1. The external sources are assembled in 
the covariant derivative, D^ = d^ + [V^, *] + {A^, *} with vector sources V^ and axial-vector 
sources A^ and in the complex density 



X = 2BoiS + tP) 



(5) 



with scalar sources S and pseudoscalar sources P. The current quark mass rriq, as usual, is 
included in the scalar source. The constant Bq is related to the quark condensate, i{qq) = 
— Bq Fq, the order parameter of DCSB in QCD and appears in the Gell-Mann-Oakes-Renner 
relation for pseudoscalar meson masses 0. 

In the large- A'^c bosonization approach, i.e. when restricting oneself to the calculation of 
diagrams with only one quark-loopQ the leading contribution to the effective lagrangian £4 
of chiral dimension 4 is parametrized [53] by nine structural constants 1^: 



£f^ 



tr 



-h (DluyDlU + h {{D^x)'D,,U + D^x{D,UY) 
+h D^U{D^U)\xU^ + Ux^) - h {Ux^Ux^ + xU^U^) 

-i,(x^u-u\)ii[x^u-u^x 

-Is (f^,d,u{d,u)^ + f;:^{d,u)^d,u) - hU^Fl^^UFl:^ 



+cr 



in Euclidean denotations. Herein Fj^^ = d^L^, — di,L^ + [L^, Ly] and 



(6) 



L 



V... 



A^] H^ 



V... 



A,. 



The constants /« {i 7^ 7) arise from the bosonization of one-loop quark diagrams in 
the soft-momentum expansion. The coefficient Ij is essentially saturated E^ by the non- 



perturbative correlator of gluon pseudoscalar densities (see |22] and references therein). The 



^This yields only terms with one connected trace operation. 



bosonization yields also vertices of dimension 4 in D™ which are invariant under chiral 
transformations of external sources, 



ninv 



-i/itr 



-pL pL _|_ rpR TpR 



-Hot! 



X^X 



(7) 
and 



The Wess-Zumino-Witten action Swzw includes the so-called anomalous vertices 0. 
it is not displayed here. 

The low-energy phenomenology of pseudoscalar mesons is described by the standard 
Gasser-Leutwyler lagrangian which contains ten structural chiral constants Li, 



L 



GL 



-Li (tr 
-Lgtr 
+1.5 tr 



\d^u)^d^u 



Loti 



{D^U)^D,u]ii\{D^U)W,U 



{D^U)^D^U{D,U)^D,U] + L4tr [{D^U)'^D^u] tr \x^U + U^ 
{D,U)^D,U{x^U + U^x)] - U (tr [x^U + f/t 



'X 



L, (tr [x^U - U^x] ) - ^8 tr [x^Ux^U + U^xU^X 



-Lgtr 

■•ir 
-4 



F;iD,u{DMy + F;:^{D,uyDM 



L^ntT 



'X 



U^F^^^UF^^" 



+cr 



We will neglect meson-loop corrections in the following^, keeping our attention to the leading 
order in l/Nc- Then the (tree-level) relations between coefficients Li and Ij are derived 
following the usual on-shell scheme of Chiral Perturbation Theory, i.e. imposing the equations 
of motion from the Weinberg lagrangian £2, 



U^Dp - (DluyU - x^U + U^x = Tptr \u^X - X^U 
where Np = 3 is the number of light flavors. One flnds, 



(9) 



2Li — L2 — h, L3 



h + h- 2/1, 



LT = h- -h + — /3, ^8 = -7/3 + -T^h + /( 
D i/ 4 / 



L^ — Lq — 0, L5 : 



h + /s, 
Lio = /g 



(10) 



Thereby the predictions of a bosonization scheme can be directly referred to physical char- 
acteristics of pseudoscalar mesons. 

The low-energy QCD bosonization yields the chiral lagrangian in the soft momentum 
expansion whereas the ChPT concept is based upon the soft momentum and light meson 
mass expansion dim[p^] = dim[m^] = 2. Therefore chiral coefficients in the bosonized action 
have different status in the canonical and chiral dimensional analysis. Namely, the canonical 
dimension of vertices containing powers of the field x (see eq. (^)) is always less than the 

^Consequently, we will not achieve higher precision than the usual logarithmic uncertainty due to the 
choice of renormalization scale. 



chiral one because the canonical dini[x] = dim[S'] = 1 and the chiral dini[x] = dini[m^] = 2 . 
The corresponding structural constants h,5fij or -^4,5,6, 7,8 and H2 are in fact dimensional in 
the canonical sense that becomes evident when to supplement them with factors of Bq from 
X- Thus these constants as well as the coefficients of £2 are essentially non-perturbative as 
they are proportional to powers of the basic QCD scale 

where g'^{fi) is a QCD gauge coupling constant measured at a scale fi. 

It happens to be that the dimensionless structural constants /i, 2,3,8,9 or -^1,2,3,9,10 are not 
sensitive to a DCSB model and take universal values [|, p!5[] , 



4/1 = -2/2 = 2/3 = /8 = 2/9 = -^. (12) 

4o7r^ 

These constants are in a satisfactory agreement with phenomenological values obtained from 
the strong and electromagnetic interactions of pseudoscalar mesons [|, |^, [l^, [19|. 

On the other hand, the values of other, dimensional chiral constants (except for 
L4 = Lg = which are suppressed in the leading l/Nc order [g]) are strongly sensitive to 
modelling the DCSB by means of a momentum cutoff A and a spectrum asymmetry M or 
a dynamical mass M^. In order to estimate dimensional constants in a model independent 
way we extend the chiral lagrangian by introducing in the next section scalar glueball, h and 
quarkonium, S variables so that their v.e.v.'s provide the required constants. Further on, 
we will normalize the corresponding lagrangian to reproduce the QCD motivated expansion 
of the vacuum energy for small and large external scalar sources. It will provide the model- 
independent constraints on coupling constants of the joint scalar-pseudoscalar lagrangian. 
Finally, the large mass reduction of scalar mesons in the latter one yields the required 
estimates on the chiral constants L^,L^. 

Let us start this derivation with preparing the concise form of the dimensional part of 

'~'chi 



CT = tr 



8Blh {S^ - P') - ^BlH2 (S^ + P')] + ACf^ , (13) 



where DY^ = 9^ + [\/^, *]. Here and in what follows we omit the Ij term as we do not discuss 
the U{1)a terms in this paper. The chiral fields U{x) are encoded in rotated external sources 
by means of the following chiral bosonization rules, 

V, + A^ = U\V^ + A^)U + U^d,U- %-A^ = V^-A^- 
S + iP = U\S + iP); S-iP={S-iP)U . (14) 



With these rules one can easily reproduce the chiral lagrangian (y). The terms in (|I3|) are 
ordered in accordance with the ChPT count and A£^™ stands for vertices of dini-6 and 
higher. Meanwhile from the canonical dimensional analysis it follows that two more vertices, 
namely 

A£^f = 85o/io S"" + SBoIn SP^ + . . . , (15) 

have the comparable dim-3 and they should be retained in the extended scalar-pseudoscalar 
lagrangian (see below). 

The bosonization models with DCSB introduced by an asymmetric cutoff A, M |Q or by 
a dynamical mass M^ [^, |I^, [T^ yield in fact more relations between the above constants. 



Let us derive them first in the framework of the chiral quark model |rT|, |T2[. We impose 
that the quark loop effective action in external fields were regularized (see [0, |28|) by 



means of a momentum cutoff A^ so that it is left (and right) invariant under local chiral 
transformations of external sources like in (0) (we ignore here the P-odd Wess-Zumino- 
Witten action) . It corresponds to the definition of chirally invariant quarks. The leading 
low energy contribution to this action in the large cutoff approximation A^ ^ l^l, \P\, \V\, \A\ 
is given by divergent parts of loop integrals |]T^, |TU|, ^ , 

r, ~ -^Jd^xtT[~2cAl{S^ + P') 



167r2 

+ ln^ {{{S' + Py-[S,P] 



+ {{Dlsr + {DlPf - {A„ Sf - {A„ Pf - 2tDlP{A„ S} + 2tDlS{A„ P} 



-l{iF^.r+iFp')] 



+ ■■■ (16) 



It consists of four independent chiral invariants corresponding to the four lines of (pig). 
For an 0(4) invariant momentum cutoff one has c = 1 but it is regularization dependent. 
Meanwhile, the logarithmically divergent coefficients are unique. 

The coupling of chirally invariant quarks to pseudoscalar fields is provided by the dynam- 
ical mass term introduced in such a way that the designed chiral invariance of (the P-even 
part of) the fermion determinant is preserved, 

{S - 275^) ^(S- tl^P + M, (u^Pl + Pr))^(S- t^,P + Md) , (17) 



following the notations of (0), and as well V^, A^ — > V^, A^. The equivalence is provided by 
the invariance under quark field rotations ql -^ Uq^. The shifted effective action ([iBD yields 
the dimensional part (|13D of the chiral lagrangian C^i^ in terms of the parameters A^^, M^, 
which simulate the DCSB, with the following predictions for the coefficients. 



^o^^S^ln^; i?oFo^^^(cA^-Miln^); ^ ^ M,- 




^ - / - ^° ■ J - 


= 0; 


h = 


Fi . 


'' ^'' 8B0M,' ^' 

^0 ^0 J 1 or 


16i?2' 


8B0M, 8BI '' ' "'° • 





(19) 

Inspecting other bosonization models IU-0] we find that the relations for Ij displayed 
in ([T9|) are also valid in the large-cutoff approximation. For instance, in the QCD chiral 
bosonization model |0] the DCSB is introduced by means of an asymmetric regulator for the 
quark determinant, 9 (A^ — (^ — iM) j which is non-invariant under local chiral rotations. 
The pseudoscalar meson fields arise from local chiral rotations of the quark fields q{x) = 
{PlU{x) + PR)qinv{x)i and the chiral lagrangian is made from the chirally non- invariant 
part, 

^= /P^ f'^'tt 5! ^^-^^^ /pt/exp(-^e//(f/;V^,A5,P))(Z,„.)G, (21) 

where {■ ■ ■)g stands for the averaging over the gluon vacuum. Seff does not contain any gluon 
fields. The role of gluons is reduced only to the formation of the dimensional parameters 
of the theory, (A, M), on the base of the equation of stability of the low-energy region, 
supported by the gluon condensate [^ (see the next section). After calculation of the chirally 
non-invariant part of the quark determinant one comes to the chiral lagrangian (^ with 
dimensional vertices collected in (p^ . The description of the main parameters and the 
structural constants in the minimal version of the model in [Q has the following form: 

F^ = ^ (A^ - AP) ; BoF,^ " ^ (^' " l^") ' (^2) 

^' = ^'' = ^'' = TdWo'' ^^ = °' ^^ = "16%' ^^^^ 

H.--§i-^h. (24) 



When comparing with (|l^) one can see that the following relations between the constants /, 



h = ho = /ii ; /5 = ; h = --^ ; (25) 

are the same in both models, but the analytic dependence on model parameters and the 
quantitative predictions for the coefficients are (of course) model-dependent. In particular, 
the values for H2 do not coincide for any choice of parameters. As to the other constants 
Fq, Bq, Ii, their values may be adjusted to the same values in both models but for an unnatural 
choice of the constituent mass parameters: Md — 100 MeV in the chiral quark model and 
M ~ 600 MeV in the chiral bosonization model [^, respectively. In the following, we will 
therefore only use the general and model-independent relations (B^) to reduce the number 



of unknown parameters in the chiral lagrangian (|T3|) or its extensions including heavy scalar 
fields 1^ (see below). 

For illustration, let us however check the consistency of these predictions with the phe- 
nomenology of pseudoscalar mesons (PUf ) for Ma — 250 MeV and Bq ~ 1.3 GeV (A^ ~ 
1 GeV): 



^5 



F'o 



3 • 10"^ 
1287r2 2 ^ 



1 



^0^ 



1652 



(cf. (2.26 ± 0.14) ■ 10"3 from [29]) ; 

0.5 • 10"^ (cf. (0.9 ± 0.4) ■ 10"^ from |T 



(26) 



where iVc = 3 is taken. Thus the chiral quark model displays a satisfactory agreement with 
experiment. Note that whereas L5 = J4 rather serves as an input parameter, which may be 
tuned by a suitable choice of model parameters, the value of Lg is well calculated from all 
the inputs. A similar situation arises in the chiral bosonization model. 



3. QCD bosonization in the scalar sector 

In the large Nc limit the main contribution of the scalar glueball sector to the effective 
quarkonium lagrangian is provided p4| by the v.e.v. of the glueball field < h >= {{G^^j^Y)n.p., 
i.e. the gluon condensate (see (H)), 



Cr. 



47]-2 



{{G%f)n.,. ^ (400 MeV)^ 



(27) 



Here n.p. means that the perturbative part has been subtracted. Let us associate the 
scalar quarkonium fields with the radial colorless fiuctuations of light quark fields, q{x) = 
exp||E(x)|- qinv{x). Herein gj„„ are fields invariant under dilatations. In general, the di- 
latations mix fiavor numbers, S(a;) = cri + ctjTj where Tj are generators of the SU{Nf) 
group. This factorization can be used in the conformal bosonization procedure to derive the 
effective scalar-pseudoscalar lagrangian [^ |2^. We evaluate it in the soft momentum region 
neglecting the kinetic terms and higher-order derivatives of (heavy) scalar meson fields. 

The dimensional part of the lagrangian for the fiavor-singlet (dimensionless) scalar field 
a then has the general form: 



C 



dim 



tr 



a 



^g^2 \(Gf^'^) /"-P- 



+ -a^e 



-4a 



+ a^S e 



-3a 



+ {a2iS^ 



022 P - 023 Af,) e 



+ (an S-' + ai2 SP' - a^ lA^D^ P - 2aiiA'^S) e"'^ 



+ ... 



(2^ 



Let us understand this form of the lagrangian in more detail: The first term simply repro- 



duces correctly the scale anomaly of the quark determinant [30|-[32|. The following terms 



are ordered according to their canonical dimension, i.e. they have dimensional vertices with 
coupling constants which are polynomials in the QCD scale Ac, CLjk ~ A^ (the dim-4 term is 
a purely scalar vertex). According to the dimensionality of the coupling constants the ver- 
tices have been dilated with appropriate powers of exp(— cr). We have again neglected dim-4 
perturbative contributions and higher dimensional vertices proportional to inverse powers of 
Ac- It will be shown in the next section that such an expansion can be made consistent with 
the basic properties of the QCD vacuum energy. 

Let us first normalize the lagrangian such that the minimum of /^^™ for vanishing external 
sources is just reached at cr = 0. It corresponds to 

"^ = 4^ ((^^^)')«-P- = § • (29) 

In order to find the chiral structural constants /4,5,6 one should perform the saddle point 
approximation for the effective scalar E-field lagrangian [^ and develop the heavy scalar 
mass and ChPT expansion. Of course, when the external sources are present, the lagrangian 
(^) does no longer possess the minimum at S = 0. The saddle point then is rather given by 

Sets ;=, , / ^21 Qttg \ 52 , 1 / 52 42 



^rmn = -r^S+[-^-—^\ S' + — [a22 P' a23A'] + ■ ■ ■ ; S = S^,„ + E. (30) 

4a4 \2a4 IGaiJ 204 ^ ^■' 

Now, in terms of the shifted variable S, the minimum occurs at S = 0. Let us therefore 
define 

Cf^iJ: = S„,) = Cf''\a = 0) , (31) 

which should be then identified with the dimensional part of the chiral lagrangian ([T3|), 
£^™(S = Sm) = C^j^. Consequently, the coefficients of the chiral lagrangian (|T3|) can be 



matched to (|28D in the following way: 



% = ^3 = -BqFq ; 

41=^21 -^ = -452 (2/6 + 1/2); ^2 = «22 = 45^(2/6 -i/2); a23 = «23 = ^o; 



504 

303021 A5al _ „ j5 . , _ 303022 

2O4 32O4 2O4 

3O3O23 



, 003021^003 n' -n ■3"3022 _ q p r 



^13 = «i3 = 8^0/4; 0^4 = oi4 = 8Bo{h + h) ■ (32) 

404 

Just these coupling constants should be compared with the chiral bosonization predictions. 
In particular, the quark model lagrangian (p!^)-([20|) (or the bosonization model ([2T| ) - 
is induced by the choice: 



"'" - ^» {' - 2I 



or = If:- 



FqBq 



^22 



"11 



''21 



''12 



''23 



''13 



2M, 



or 



-H 
2^°. 



''14 



F2 



or 



27r2 



(33) 



In the next section we will remind basic properties of the QCD quark vacuum energy as a 
function of external scalar fields and further derive the additional constraints on some of the 
ttjk- It will allow us to reduce the number of input parameters in any model building. 



4. QCD vacuum energy and large scalar fields 

We go back to the QCD generating functional (|l|) and consider nearly constant scalar gluo- 
nium p and quarkonium S sources. The renormalized coupling constant at a given scale n 
for Np active quarks obeys the following RG equation 

1 , , /i 11Nc-2Nf 



drg' =(3{g' 



-bo 9 ; 



g'if^) 



bn In 



Af 



247r2 



In 



An 



(34) 



in the 1-loop approximation. Herein dr = (9/(9 In (/i/ Ac). 

From eq. ([^) it is obvious, that the field p can be reabsorbed into the definition of the 
coupling constant, i.e. the definition of the basic QCD scale Ac, 



11 

:r^ = ^ + P = Oo In ■ 



/^ 



-p/bo 



Ac = Ace" 



(35) 



g- r Ace- 

Now let us perform an expansion of the vacuum energy in Np, the number of (light) flavors, 

£Up, s) = e^^lip) + Np e'-^lip, s) + 0{nI) . (36) 

The leading term for Np = is purely gluonic. The non-perturbative part of the QCD 
gluonic vacuum energy density £^^^{p) is then defined as 

with fi = / d'^x being the space-time volume. 
From dimensional arguments we obtain that 

1 



(37) 



£liP = 0) = (644) 



n.p. 



■(©M.) 



Hfi/n.p. 



-A'c, 



(38) 



where 0^,^ stands for the energy-momentum tensor. One therefore arrives at the RG equation 
for the energy density dr S (0) = —4 6 (0) which has the following (RG invariant) solution 
il|] (see eqs.(|l|), (|3^) ) given by the trace anomaly [p0|-p2|: 



^lac(o) = Y^/5(/)((G;:jVp. 



Np=0 



16 



{(G^^u) )n.p. , 



(39) 



10 



where bo = 60 I at =o- When p 7^ one simply has to replace Ac with Ac from (|35D 



^liP) = ^-'"~'' £l?c(0) = -^ C, e-^'/'^o , (40) 

in terms of (|27D . 

In order to evaluate the effective potential for the scalar glueball field h{x) we follow 
the bosonization ansatz (^ which can be rewritten in the form of a functional Legendre 
transform with respect to the (nearly constant) field p, 



Z{h) =exp{-fiyg(/i)} = jVp exp|-fi 
For f2 — >• 00 the saddle point configuration, 



'e^'Hp)--~ph 



(41) 



'=^«''^-T'"Wk;' <''' 



delivers the required effective potential: 



h , /i 



^.(^) = <ai(0) TT 1 - 1^ TT ' ^0 = ((G';:.)')n.p. , (43) 



which coincides with the expression in |p2| . 

Now we proceed to the analysis of the quark- loop contribution in (|l|) to the vacuum energy 
and its dependence on the external scalar source S{x) ~ const (other external sources are 
not shown for the time being), 

Z(p, S) = exp {-Vt e^ac(P> S)] 



JvGe^vl^-J d'x ^ {Gl^f + AT^ 1 Tr [in ^E±ll±^']^ r 



> , (44) 

n.p. 



where the quark fields are integrated out leading to the quark determinant. Here R = 
R{p'^,A'ljy) is a gauge invariant regulator. In the following we restrict ourselves to slowly 
varying sources ^5* ^^p ~ 0. 

The first-order term £ gets two contributions PH] from both the /^-function. 



,2 



5..r^--^ln^, (45) 



1 , /ij 

and the quark determinant. In addition, all dimensional operators are supplemented with 

scahng factors exp(— p/60) so that the higher- dimensional condensates^ C2„+4 ~ 

^We retain the 1-loop contributions only and therefore neglect the anomalous dimension of qq{x), 
i.e. of S{x). 
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(Tr [^2n^]) 



n.p. 



Af""*" are rescaled according to their canonical dim = 2n + 4 and (|35|). 



Thereby one derives that 



^l.>|(Gy'-itr 



{x\ In 



>2e-2p/feo + ^2^ 



/i^ 



i?|x) 



J / n.p. 



(46) 



It is a RG invariant quantity [p5l , 



9£ 



(1) 



vac p— 4p/bo 



5 In/i^ 
which is provided by the Fujikawa's theorem |3^ 



^(G;:.)^ + ^tr[(a:|i?|x)] 



n.p. 



{tr[{x\R\x)])^^,^ 



\\^fiu) /n.p. ; 



(47) 



'n.p. 48^2 

as long as the regulator obeys the conditions R{P^ = 0) = 1, R{P^ — > oo) = (in fact any 
momentum cut-off scheme). 



Let us introduce a fermionic reference scale A^ [^ for the Dirac operator P as a nor- 
malization at S* = in the following way: 



tr 



(x\ In --7^ R\x) 



= 



leading to 



£t'MS = o) 



3-4p/feo 



1 

24 



G 



in.^^ 



2p' 



Ah bo 



C„ 



(48) 



(49) 



Presumably, this (fermionic) scale is of the order of 1 GeV, Kp > Ac- Then we can redefine 
the quark vacuum energy using this scale, 



e^^^{p,s) = £'-^^{p,o) + e-'^rb^Wp{S) 



WpiS) 



-- ( tr 
2 \ 



{x\ In — -2 R\x) 



A2 



(50) 



n.p. 



where we have defined S = 5*6^/*°. When other external sources are present, they should be 
rescaled in the same way, according to their canonical dimension, (P, D^, A^) = 



The glueball potential is again given by the Legendre transformation (HTI), and to the 



first order in Np the quark contribution is described by (50) where for the field p the saddle 
point value p := p{h) ( ^2] ) is taken. These rules can be straightforwardly applied to QCD 
bosonization models in order to derive the joint scalar glueball-quarkonium potential. 



exp 



~nwF{S)e- 



-ip{h)/bo 



T>a exp 



.Q [C^^S)-^^]e~'^^'^^/'° 



(51) 
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The constant in the exponential on the r.h.s. is fixed by the normahzation Wf{S = 0) = 
provided at the saddle point a = 0, and for our purposes £cr is taken from eq. (|28|) . After 
substituting (^Tj) into the QCD generating functional (H) and performing the shift a -^ 
o" — p/bo, the glueball-quarkonium potential may be eventually expressed in the following 
form: 



A.^^.(/.) + ^/.(a-ln^ 4 



\) + ^f^iS) 



12 



a 



(52) 



We see that in (^2|) the coupling between glueball and quarkonium fields is fixed unambigu- 
ously. 

Let us examine the dependence of the QCD vacuum energy on slowly varying scalar 
fields, S{x) ~ const. The quark vacuum energy is connected to the quark condensate (at 
1-loop approach) [|^ 



dsWF{S) =i{^>m) 



tr 



{x\ 



,2—^R\^) 



1Z)' + S^ 



n.p. 



sign(S)(tr {x\nS{p^) R\x) )„ 



S 



(tr[(a;|i?|x)])„.p. 



US^' 



(S^oo) 



(53) 
(54) 



Let us focus our attention on the case S ^ Ac- From eq. ( PI) we have 

Z(p,5^oo)^|PGexp|-|d^a;i(^l-iV^^ln^^ (G%Y + OiN^f 



n.p. 



(55) 



As in 



one can make the RG improvement by redefining the strong coupling constant 

xxiv ju^ o/\r_ C2 1 1 /\r ,,2 

- - ^TT^ In - 
/i 



2 ~-2 

9 =9 



Nt? In — ~ 



UN, /i 

T In ^^— 

487r2 A2 



c 



2Nf , S 

t: In ^^— 

487r2 A2 



c 



lliV, ^r_ 

487r2 A2 ' 



(56) 



so that 



Ac = Ar 



.An 



(57) 



Consequently, the leading contribution into vacuum energy for large S factorizes. 



-(0), 



'\S\' 



SNp 
lliVc 



£v^{S -> oo,p) = £;^^{p) \j- 



+ 0{Nl) < . 



(58) 



In all cases quark polarization effects lower the QCD vacuum energy. 

Passing on we remark that this result is related to conclusions in ref. |33]. Indeed let us 
assume the value of the strong coupling constant and therefore of Ac is fixed at a scale much 
larger than the masses of all Np = 6 participating quarks (say e.g. a GUT scale Mu). Let 
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us denote this scale by Ac- If one likes to know how this is connected to the low-energy Ac 
with Np = 3 dynamical quarks, one has to integrate out the heavy quarks step by step. In 
this way one extracts the large mass logarithms in a similar manner. The result is 

/ \ 2/27 

Ac = Ac f !^^£^ j > Ac; 27 = IIN^ - 2N^''' . (59) 

Thus in this Sec. the structure of bosonized action for a scalar glueball and its coupling to 
quark degrees of freedom is determined as well as the asymptotic conditions on the behavior 
of quark vacuum energy for very large but constant scalar sources 5* are formulated. We 
would like to stress that this limit belong to the soft momentum region and it is safe being 
related to the decoupling of heavy quarks, i.e. one does not expect any phase transition or 
strong coupling phenomena when S* — > oo. Hence one should use this conditions to obtain 
model independent constraints on parameters of effective meson lagrangian. 

5. Asymptotic constraints on chiral constants 



Let us impose the QCD asymptotics (p5D to be fulfilled also in the extended meson lagrangian 
(PSP for scalar and pseudoscalar fields. We normalize a^ according to (PD| ) and 03 as in (P^). 
First we switch off all the external sources but the scalar one, i.e. we retain the scalar field 
vertices with constants 021,011. 

For these unknown parameters, we examine the asymptotics 5 ^ 00 {V,A,P = 0, 
U = 1). It is convenient to introduce a new variable y such that 

Sy:=e-'' ^a = -\nS-\ny . (60) 

We then obtain for the effective potential 

W{y,S) = -04 InS"- a^lny + {-a^y^ + a^y^ + 021!/^ + cmy) S'^ . (61) 

The first term (—04 In 5) already reproduces the leading logarithmic behaviour (^) of the 
QCD vacuum energy with the correct coefficient. Consequently, after integrating out the 
field y from the saddle point approximation, the remaining terms should at most behave like 
(const. + 0(1/^2)). 

Thus the term in front of S'^ has to be of higher order in l/S"^ for y taking its saddle 
point value y^m 

ymin f ^04 yl,i^ + 03 yl,i^ + 021 ymin + ^n j =0 ( — ) ] ymin ¥" • (62) 

For ymin being the saddle point, the first derivative of W{y, S) with respect to y has to 
vanish, 

-I 

W{y) = a^ yfriin + 3 03 y^^^ + 2 021 y^in + oii - ^ = . (63) 
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Let us expand the last relation in 1/S'^, with ymm = Vo + S/S'^ + 0(1/5"^). We obtain 



W'{y) = 04 yo^ + 3 as 1/0^ + 2 021 z/o + «i 



5 (3 04 1/0^ + 603 2/0 + 2a2i) , ^. 1 ^ .„,s 

+ ^1 + 0{—). (64) 

Hence one arrives to the second constraint, 

04 yo^ + 3 as yo^ + 2 021 ?/o + an = , (65) 

Two constraints make it possible to find the saddle point value yo and to estimate one of the 
coefficients, an: 



— 4as ± Y 1603 — 1202104 

yo = 5 ; 

004 
8a| - 604021 , 4os02i ,„„. 

""'' = ^° 9Z + ~9^- ^^^^ 

In order to obtain information for the chiral lagrangian, namely the constants L^jLg, we 
perform the large scalar mass reduction Oj(j) -^ K(j)y see eq. (|32D . and express the chiral 
constants in terms of Bo,FQ,Cg and Og^. Together with the universal relations (|33|) a[i = 



a'i2 = fl'is = 014 = 8-00-^4, h = one finds the parametrization of the chiral constant L5 






9 L 



4 16 V2 



^-^- («^) 



The chiral constant Lg is supposed to obey ( PU]) and will be in agreement with the exper- 
imental estimates for the choice of L5 ~ 3 ■ 10~^. Note that the pre-factor in (|67D already 
sets the right scale {BqFq /?)C'^] ~ 10^'^. Therefore the remaining numerical factor in paren- 
theses is expected to be of order 3, which is in fact only achieved for the choice L'^ . It 
can be achieved for a value of O21 ~ ^\Fq < which just lie inbetween the predictions of 
the chiral quark model a'21 = Fq (1 — {Bq/2Mii)) ~ —(2 -^ 3)Fq and the chiral bosonization 
model O21 = 3/2 Fq, respectively. Thus we see that the asymptotic constraints impose rigid 
conditions on the choice of a bosonization model. 

The next-to-leading term of the asymptotics (|6^) needs a special care, li 6 = then the 
next term of order 1/S^ contributes to the quark vacuum energy ( |6lD only corrections of the 
same order in virtue of (|63D. Then 0(1/ S"^) terms in the effective potential do not appear. 
That would correspond to the identically vanishing dim-6 gluon condensate |^ Cq ~ (G^)- 



As there is no reason a priori to neglect the latter one we obtain the third constraint for 
MO, 

3o4i/o^ + 6031/0 + 2o2i = . (68) 
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Together with (^), ( |65D we thereby have 3 constraints for the variables a2i,cin,yo- The 
analysis of these 3 constraints shows that their compatibility holds only if 021, an 7^ 0. One 
then has, 

1 « D 172 1 « d2 174 op;«d3i76 

(69) 



IQBoF^ IQBiFi 
Vo ^ J «2i r^ ' '^11 


25653^6 


which corresponds to 




^^ = ^t = L^ = §^- O(IO-) ; 


-1- 



(70) 



However, the amount of terms in the polynomial expansion of (|28|) may still not be enough 
to make a numerically reliable interpolation, and it will be necessary to include terms of 
0{S^) into the scalar lagrangian (^81) . (Note that the unnaturally large value of an which is 
the last coefficient in the polynomial expansion of ( pSj ) will presumably reduce then, leading 
also to more realistic values for L5.) Therefore, the third constraint ( |68D cannot be taken too 
seriously for a given interpolation, and neglecting the dim-6 condensate corresponding to the 
solution 6 = may well be consistent with the precision of the expansion of the bosonization 
model (|2|), 

Nevertheless, we have seen that the asymptotic constraints can be embedded in the 
meson sclalar-pseudoscalar lagrangian, giving a powerful tool to estimate some of the chiral 
coefficients. 

Evidently the constants in other vertices of the lagrangian (^) can be treated in the 
same manner. For instance, one may expect the mutual cancellation of large S contributions 
in certain combinations of vertices a22, ai2 or a23, ai4, 

a22ymin + ai2 ^ 0; a23ymin + 2 ai4 ~ , (71) 

which however is not anticipated to be accurate as only two terms of corresponding expan- 
sions are included. 

6. Summary and conclusions 

In this paper we have developed the modified bosonization approach to construct meson 
lagrangians from QCD where the large field but low-energy asymptotics in QCD is embedded 
after bosonization for scalar field vertices. 

First, it was shown that the chiral bosonization in different approaches results in the same 
structure of vertices (in the large-cutoff approximation), encoding the information about the 
QCD dynamics in a few chiral constants Jj which can be related to the phenomenological 
structure constants Li and obey certain model-independent relations among each other. 
Next, we have determined the vacuum properties of the QCD quark energy in the presence of 
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large external scalar fields and applied the results to estimate some of the coupling constants 
in the extended effective meson lagrangian, including scalar and pseudoscalar fields. These 
ones in turn have been used to obtain an unambigious estimate of the phenomenologically 
interesting (but generally rather model-dependent) structure constants L5 and Lg. 

Moreover, by employing conformal bosonization methods one gets a scalar meson la- 
grangian which satisfies all QCD motivated requirements and asymptotics and thereby gives 
a direct way to extend it unambigously, such that the coupling to glueball fields can be 
described. As an interesting result, we have given a path-integral derivation of the effective 
glueball potential from the generating functional of QCD using (5-functional constraints. 

There are a few improvements to be done in order to reach better precision, and one of 
the most important is to include higher- dimensional terms into the polynomial expansion 



of the lagrangian (pS]). In a sense this will lead to an infinite number of asymptotic sum 
rules for an infinite number of unknown constants, relating them to gluon condensates of 
increasing dimension. 

Finally, note that in the conformal bosonization approach dynamical scalar mesons are 
treated as intrinsic dilatational modes. It is a challenge to introduce them in the QCD 
effective action as true collective variables, in analogy to chiral rotation modes describing 
pseudoscalar mesons. This program will be realized elsewhere. 
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